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ABSTRACT 

We present optical integral field unit (IFU) observations (VLT/VIMOS-IFU and 
WIYN/SparsePak), and associated archival deep Ha imaging (MPG/ESO 2.2m WFI), 
of the nearby starburst galaxy NGC 253. With VIMOS we observed the nuclear 
region and southern superwind outflow in detail with five pointings, and with the 
WIYN/SparsePak IFU we observed two partially overlapping regions covering the 
central disk and northern halo. The high signal-to-noise of the data and spectral res- 
olution (80-90 kms -1 ) enable us to accurately decompose the emission line profiles 
into multiple components. 

The combination of these datasets, together with the wealth of information on 
NGC 253 available in the literature, allow us to study the starburst-driven superwind 
in great detail. We investigate the known minor axis outflow cone, which is well-defined 
in the Ha imaging and kinematics between radii of 280-660 pc from the nucleus. 
Kinematic modelling indicates a wide opening angle (~60°), an inclination consistent 
with that of the disk, and deprojected outflow speeds of a few 100 kms -1 that increase 
with distance above the plane. The [Nil] /Ha and [Sii]/Ha line ratio maps imply 
that a significant fraction of the wind optical emission lines arise from shocked gas, 
with localised pockets/filaments of strongly shocked gas. From the kinematics, the 
cone appears partially closed in at least one place, and very broad Ha line widths 
(>400 kms -1 FWHM) suggest there is material filling the cone in some regions. 
Extrapolation of the cone to its apex shows it is not centred on the starburst nucleus, 
suggesting the wind is deflected and collimated by the dense circumnuclear material. 
We discuss the implications of these findings on our understanding of the origins and 
evolution of the superwind. No evidence for an outflow is found on the north-western 
side of the disk out to >2 kpc in our optical data, due to obscuration by the foreground 
disk. The lack of an obvious connection between the inner (r<l kpc) Ha and X-ray 
bright outflow cone and the large-scale (r<10 kpc) X-ray "horns" is also discussed. 

Key words: galaxies: individual (NGC 253) - galaxies: starburst - galaxies: ISM - 
ISM: kinematics and dynamics - ISM: jets and outflows. 



1 INTRODUCTION 

The consequences of a starburst can dramatically affect the 
evolution of the host galaxy through photoionization and 
mechanical/chemical feedback, and (if powerful enough) the 
intergalactic medium through the escape of ionizing radia- 



* Based on observations collected at the European Organisation 
for Astronomical Research in the Southern Hemisphere, Chile un- 
der programmes 060.A-9123 and 081.C-0808. 
f E-mail: msw@star.ucl.ac.uk 



tion and the products of nucleosynthesis. By studying local 
examples of starburst-driven outflows, it is possible to inves- 
tigate the underlying physics of the feedback processes, and 
then apply this knowledge to gain a better understanding 
of the phenomenon of starburst outflows. This is important 
since it is now recognised that starburst-driven outflows are 
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ubiquitous at high redshift (Pettini et al. 2001 Conselice 



Martin & Bouche 



2009). In re- 



cent years, observations of local systems have been providing 
more and more evidence that superwinds are not the sim- 
ple, smooth, thermally-driven outflows that were once envis- 
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Figure 1. MPG/ESO 2.2m WFI continuum subtracted Ha: image of NGC 253 (square-root scaled and inverted for clarity) with the 
VIMOS-IFU (green boxes) and SparsePak footprints (arrays of red circles) shown. Dashed lines indicate the extent of the projected x\ 
and X2 bar orbits ( |Sorai et al.| 2000 Das et al. 2001). Inset shows a zoomed out version of the same Ha image with a box outlining the 
field-of-view of the main figure. 



aged (e.g. |Chevalier & Clegg|1985| ). A starburst is typically dra ( |Strickland k Stevens] |2000| ). |Strickland et al.| (2002) 



very clumped, with each clump containing multiple super- 
star clusters (SSCs). This seems naturally to lead to highly 
structured, multi-phase (10 2 -10 7 K), mass loaded winds. It 
is not known, however, how the flow from individual star- 
cluster sources is collimated, where it is accelerated and im- 
portantly if /where mass is entrained/ loaded into the flow. 
To address these issues and others, we have studied the inner 
wind regions of the nearby starbursting galaxies NGC 1569 
and M82 using optical integral field unit (IFU) observations 
( Westmoquette et al.||2QQ7a|b[ [2Q09a|b| ). 



In this paper we extend this project to examine 
NGC 253, a much larger system with a nuclear starburst 
and large-scale superwind. We present optical IFU obser- 
vations and associated deep Ha imaging in an attempt to 
understand the inner region of the outflow in more detail. 

Early, deep optical and X-ray imaging of NGC 253 
showed the presence of an extended (r < 1 kpc) nuclear 



outflow (Demoulin & Burbidge|1970 


Ulrich|1978 Fabbiano 


& Trinchieri||1984 McCarthy et al. 


1987). This superwind 


has been observed in optical emission lines ( Schulz &; Wegner 



found a strong spatial correspondence between the {Chan- 
dra) soft X-ray emission and Ha emission on scales down 
to <20 pc, implying that the standard adiabatic model of 
superwind flows (Chevalier & Clegg 1985) cannot hold. In- 
stead of tracing the hot component of the wind outflow as 
the model predicts, the X-ray emission, like Ha, primarily 
traces the interface zones between the hot and warm phases 
of the outflow. 



In addition to this nuclear outflow, huge lobes of dif- 
fuse X-ray and far-IR dust emission extend up to 8-9 kpc 
away from the disk, and are thought to trace the walls 
of two large-scale bubbles containing low-density hot gas 
dPietsch et al.||2000| iStrickland fe Stevens||2000| JStrickland 



et al.|2002[|Bauer et al.|2008||Kaneda et aL|2009| ). The link, 
however, between the inner (conical, limb-brightened) wind 
studied in this paper and seen with Chandra, and the large- 
scale (^10 kpc) diffuse bubbles/lobes seen with the ROSAT, 
XMM-Newton and AKARI telescopes, is still unclear. The 
inner superwind evidently is driven by a burst of star for- 
mation that has been taking place in the central ^500 pc 
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Figure 2. (a) Zoom-in of the WFI Ha image (again greyscale inverted and stretched to enhance low surface-brightness features) with the 
VIMOS-IFU footprints shown. The dashed lines enclose quadrant 3 of the detector for which instrument problems resulted in the blurring 
of the spectral line profiles, while the white ellipse marks the nuclear region including the radio, X-ray and IR peaks; (b) HST F814W 
(/-band) image of the central region of NGC 253 (again in inverted greyscale) showing the VIMOS-IFU footprints and the SiO ((v=0, 
J=2— >1) integrated intensity contours from |Garcfa-Burillo et al.]|2 000) that trace the inner and outer molecular tori in the circumnuclear 
disk; (c) Chandra soft X-ray image (inverted greyscale) from Strickland & Stevens ( 2000 ) again with the VIMOS-IFU footprints overlaid 
(the dashed curved lines delineate the region in which [Strickland fe Ste vens ( 2000]) identify and model the limb-brightened cone outflow). 



for th e past 20-30 Myr flRieke et al.119 80 ; Engelbr acht et al. 
1998| at a rate of -3 M Q yr" 1 < |Ott et al.|2005| ~ 



The nuclear starburst region itself is highly obscured 
(Ay ~ 5-18 mag; [Kornei & ; McCrady||2009| |Fernandez- 



Ontiveros et al. 2009) and morphologically complex. The 



resulting confusion has caused a proliferation of conflicting 
views of the nature of the nuclear region. On the smallest 
scales, optical, near-IR and mm sources appear to trace out 
a ring of size -50 pc ( |Forbes et al.|2000| ). This region is con- 
tained within two nested circumnuclear tori observed in CO, 
SiO a nd H 13 CO + | Garcia- Burillo et aI.|2000||Paglione et al. 



|2004| ), and this, in turn, within a much larger (—500 pc) 
cold gas torus associated with the 7 kpc long stellar bar 



(Scoville et al. 



1985 



Mauersberger et al. 



1996 



Das et al. 



2001). The galactic nucleus of NGC 253 has historically 



been associated with a powerful, compact, non-thermal ra- 
dio source (TH2) often cited to be a low- luminosity AGN 
(Ulvestad & Antonuccil 119971 I Weaver et al. 1120021). How- 



ever, high-resolution radio (Brunthaler et al.|2009|, near-IR 



( Fernandez- Ontiver os et al. 2009), and IR (Engelbracht et al. 



1998) imaging do not show any evidence of a suitable point 



source that would indicate nuclear activity. Recent measure- 
ments of the 2D stellar kinematics with VLT/SINFONI by 



Miiller- Sanchez et al. (2010) identify the kinematic centre to 



lie within a — l'/2 region containing both TH2 and the central 
point-like X-ray source (X-l). Water maser emission has also 
in the past been associated with the presumed nuclear radio 
source TH2 flNakai et aLp 995). However, [Brunthaler etaE] 
( 2009 ) showed that the 22 GHz water masers are not related 



to AGN activity, but almost certainly associated with star 
formation. Thus, although some evidence for central black 
hole activity remains ( |Turner fc Ho|1985| [Mohan et al. 2002 



Weaver et al.|2002| ), the overwhelming consensus points to- 



wards a purely starburst origin for the radio-X-ray emission 
seen in NGC 253. 

A number of young (20-30 Myr) massive star clus- 
ters have been discovered within the inner star-forming ring 
Watson et al. 1996). More recently Kornei & McCrady 



2009| ) presented IR observations of the near-IR peak, 4 
from TH2 on the south-west side of the ring, and find a 
complex system, presumably consisting of several epochs of 
star formation as evidenced by the simultaneous presence 
of RSG, OB and WR star signatures. By comparison to 
population synthesis models they find a very young age of 
^6 Myr, a mass of ~10 7 M , and Ay ~ 18 mag for this 
region. Fernandez- Ontiveros et aT] ( |2009| detect multiple 
pc-sized knots in their VLT/NACO adaptive optics imag- 
ing, and conclude that they are probably young (<10 Myr), 
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Figure 3. The extent of the large-scale X-ray superwind seen 
with ROSAT i s shown in contours jPietsch et~aTl|2000| |Strick-| 
|land et al.|2002] > overlaid on the WFI Ha; image. The dashed box 
around the nuclear regions indicates the coverage of Fig. [2J1 and 
also the approximate extent of the inner superwind. 



dust-embedded (Ay ~ 5-11) star clusters similar to those 
found in NGC 5253 and NGC 4038/9. 

NGC 253 also shows some large-scale peculiarities. 
Mal in &; Hadley| ( |1997| ) found an extended but relatively 
smooth, low surface brightness, stellar halo with an asym- 
metric extension to the south which they interpreted as 
evidence for an ancient merger. The Hi 21 cm line study 
by |Boomsma et al.| ( |2005| ) revealed two further anomalies: 
an Hi disk that is considerably smaller then the flattened 
stellar envelope found by Malin & Hadley, and extrapla- 
nar Hi around the eastern half of the galaxy extending to 
vertical distances of approximately 12 kpc above the thin 
star- forming disk and accounting for 3% of the total Hi 
mass. They interpret the extraplanar Hi as a product of 
the starburst which has lifted interstellar material out of 
the disk. Evidence that interstellar gas is being entrained 
by the wind and lifted into the NGC 253 halo also comes 
from the AKARI detection of thermal dust emission from 
the northern wind boundaries ( Kaneda et al.|2009 ). 

The structure and dynamics of the NGC 253 disk and 
halo are clearly complex. They vary both on the large and 
small scales, and are complicated by the inflows and out- 
flows caused by the bar and starburst, and by the magnetic 
field (Heesen et al. 2009a b). Unlike in M82's case, the op- 
tical outflow in NGC 253 has received comparatively little 
attention despite its prominence and obvious importance. 
We have, therefore, obtained a set of high resolution, inte- 



gral field, optical spectra of the inner regions of NGC 253 at 
a range of spatial scales, but comparable spectral resolution. 
With the VLT/VIMOS-IFU we observe the nuclear region 
and southern outflow in detail with five pointings, and with 
the WIYN/SparsePak IFU we observed two partially over- 
lapping regions covering the central disk and northern halo. 
These two datasets together provide a wealth of information 
on the state of the ionized gas, and present a unique oppor- 
tunity to study the T ~ 10 4 K phase of the inner wind and 
disk-halo interface. 

In Section [2] we present the three sets of observations we 
have used to study the NGC 253 outflow: deep Ha imaging 
with MPG/ESO 2.2m WFI, and optical IFU spectroscopy 
with VLT/VIMOS and WIYN/SparsePak. In this section we 
describe the observational details and data reduction steps 
undertaken. Section [3] contains a description of the results 
obtained from these observations, including the emission line 
maps and derived measurements. In section [4] we discuss 
what these results show, and present a simple kinematic 
model of the inner wind cone. A summary of our findings is 
given in Section [5] 

We adopt a distance to NGC 235 of 3.9 MpcQ mean- 
ing 1" ~ 20 pc, and a heliocentric systemic velocity of 
250 kms~ 1 ( |deVaucouleurs et al.|l991 ). 



2 OBSERVATIONS 

2.1 MPG/ESO 2.2m WFI imaging 

Deep, narrow-band Ha(+[Nll]) images, together with corre- 
sponding narrow-band continuum images, of NGC 253 ob- 
tained with the MPG/ESO 2.2m WFI instrument ( |Baade 
et al. 1999) were retrieved from the ESO archive to com- 



plement our IFU spectroscopy (described below). These 
were observed on the night of 21st September 2001 as part 
of programme 60.A-9123, and consist of five exposures in 
each filter, dithered to cover the gaps between the eight 
CCDs of the WFI instrument. The exposure times were 
5x900 sees in Ha (NBHALPHA/7_ES0856) and 5x450 sees 
in the corresponding 'off-band' filter (NB665/12_ES0858). 
The ^75 A FWHM of the Ra filter means that emission 
from [Nn]AA6549,6583 is also included. 

Images were processed with the imaging pipeline pack- 
age thelJ^] the steps included bias subtraction, flat field- 
ing, weighting (to identify chip overlaps and cosmetic de- 
fects), determining the astrometric solution via catalogue 
matching using SCAMP, background matching, and finally 
co-addition of the individual exposures. The on- and off- 
band images were then aligned using their World Coordinate 
System (WCS) solutions and the off-band was subtracted to 
leave a pure emission-line Ha(+[Nn]) image. Fig. [I] shows 
the full image in the inset, and the central disk in the main 

1 We note that discrepancies exist in the literature with regards 
to the distance to NGC 253; some authors (e.g. |Alonso-Herrero| 
|et al.||2003| |Engelbracht et al.||1998|> adopt the value o f 2.5 ± 
0.4 Mpc determine d by|de Vaucoule urs ( 1978 ), others (e.g. |Kornei| 
|fe McCrad y 2009 F ernandez- Ontiveros et al ||2009| ) adopt 3.9 ± 
0.4 Mpc measured by Karac hentsev et al.| (|2003| , while others 
(e.g. |Sakamoto et al.|2006| jBrunthaler et al.|2009| [Matsubayashi| 



et al.|2009fc adopt 3.5 ± 0.2 Mpc from 



Rekola et al.| (12005 1 
http://www.astro.uni-bonn.de/~mischa/theii.htnil 
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figure. Fig. [2^l shows a portion of the image covering the 
south-western outflow with a greyscale stretch designed to 
enhance the low surface-brightness wind features. 



2.2 VLT/VIMOS-IFU spectroscopy 

Queue-mode observations of NGC 253 were made using the 



Table 1. VLT/VIMOS-IFU observation log 



Visible MultiObject Spectrograph (VIMOS; Le Fevre et al. 



2003 ) Integral Field Unit (IFU) between August and Decem- 
ber 2008, covering the nucleus and minor axis outflow with 
five positions at a position angle of 140°, with 0.6-0.9 arcsec 
seeing (11-17 pc at the distance of NGC 253). 

The IFU positions were chosen to cover the nuclear re- 
gions of NGC 253, and the known inner minor axis wind 
outflow. Fig. [I] shows the position of the IFU fields on the 
MPG/ESO 2.2m WFI continuum subtracted Ra image of 
the whole disk. The location of the IFU positions are shown 
m more detail in Fig. [2] on (a) a zoom-in of the WFI Ha 
image; (b) an archive HST / A CS F814W image; and (c) a 
Chandra soft X-ray image from Strickland & Stevens ( 2000 ) . 
Fig. [2J3 also shows the SiO (v=0, J =2—^1) integrated in- 
tensity contours from | Garcia- Burillo et ah] ( |2QQQ[ see Sec- 
tion |4~2). Fig. [3] shows the extent of the large-scale X-ray 



outflow as s een with ROSAT ( |Pietsch et al.| 2 000 [[Strickland 



et al.|2002 ) overlaid on the WFI Ra image (for a discussion 



of the connection between the inner wind and and X-ray 



lobes see Section 4.5) 



Positions 1 and 2 were observed in the higher resolution 
; /33 spaxel -1 setting giving a field-of-view (FoV) of 13 x 
13 arcsecs (^250 x 250 pc). The remaining positions were 
observed in the 0'!§7 spaxel -1 setting, giving a FoV of 27 x 
27 arcsecs (-500 x 500 pc). 

The use of the high-resolution (HR) red grism (for all 
positions) gave a spectral coverage of 645-860 nm at a reso- 
lution of R=3100. This allowed us to cover the nebular diag- 
nostic lines of Ha, [N 11] AA6548, 6583 and [S 11] AA6716, 67310 
With the medium- and high-resolution grisms on VIMOS, 
the spectra are long enough that the detector no longer has 
enough space to accommodate the full 6400 spectra result- 
ing from the 80 x 80 fibres of the IFU. The IFU head is 
therefore partially masked by a shutter, leaving a square of 
40 x 40 fibres giving 1600 individual spectra. These spec- 
tra are reformatted by the optical fibres for imaging onto 
four 2kx4k EEV CCDs. These CCDs are arranged to form 
four quadrants of a rectangle and are hereafter referred to 
as Q1-Q4. 

We took a set of dithered exposures at each position 
(offset by integer spaxel shifts) . The exposures for positions 
3, 4 and 5 were split into 3 subsets of exposures. Total in- 
tegration times varied between 2340 and 6993 sees. Table [1] 
lists the dates, coordinates and exposure times for each sub- 
set of exposures for each position. 

In order to remove the pixel-to-pixel sensitivity differ- 
ences, and enable the wavelength and flux calibration of the 
data, a number of bias frames, screen-flats, twilight flats, 
arc calibration frames, offset sky exposures (exposure times 



3 The Can triplet at A8498.0, 8542.1, 8662.1 were also present in 
our wavelength range, but this part of the spectrum was affected 
by fringing which we could not correct for. 
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2x200 sees for each position), and observations of photo- 
metric standard stars, were taken contemporaneously with 
the science fields. 



2.2.1 Reduction 

Basic data reduction was performed following the standard 
ESO reduction pipeline tasks implemented in GASGANO. 
Briefly the steps go as follows: a master bias frame is cre- 
ated using the task vmbias, then vmifucalib is run to (1) 
create a master flat field frame; (2) trace the position of 
each spectrum on the CCD from the screen- flat; (3) de- 
termine the fibre-to- fibre transmission functions; and (4) 
determine wavelength calibration solutions. The standard 
star and science observations are then processed using vmi- 
fustandard and vmifuscience, respectively, to apply the 
aforementioned correction functions and extract the indi- 
vidual spectra to produce four data files per exposure (sci- 
ence and sky), each containing 400 reduced, flux- calibrated 
spectra corresponding to one of the four quadrants of the 
detector. 

The final steps of the reduction procedure were to clean 
the cosmic-rays using L ACOSMIC ( van Dokkum|200l| , mosaic 
the data files for each quadrant into one, scale and subtract 
the sky frames from the corresponding science frames, cor- 
rect each cube for the effects of differential atmospheric cor- 
rection (using an iRAF-based code written by J. R. Walsh 
based on an algorithm described in Walsh &; Roy| 1990), 
and combine the three dithered (offset) exposures for each 
position. These steps were achieved using a combination of 
iRAF0and idl. 

Due to an intermittent problem (apparently now re- 
solved) with VIMOS that caused the focussing of light onto 
Q3 of the detector to fail at particular rotator angles (un- 
fortunately including the one used in our observations), the 
spectra in Q3 of every observation were blurred and there- 
fore unusable for the analysis of line profile shapes. However, 
we retained Q3 throughout the reduction procedure so that 
we could measure the integrated line fluxes, which remain 
unaffected by the focussing problems. 

In order to determine an accurate measurement of the 



4 The Image Reduction and Analysis Facility (iraf) is dis- 
tributed by the National Optical Astronomy Observatories which 
is operated by the Association of Universities for Research in 
Astronomy, Inc. under cooperative agreement with the National 
Science Foundation. 
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instrumental contribution to the line broadening, we fitted 
single Gaussians to isolated arc lines on a wavelength cali- 
brated arc exposure in all 1200 properly focussed apertures. 
The average instrumental broadening (velocity resolution; 
FWHM) of the final dataset (not including Q3) is 2.0±0.1 A, 
or 90±0.5 km s _1 , consistent with the documented grism res- 
olution. 

The emission line profiles of Ha, [Nil] and [Sll] were 
fitted with multiple Gaussians as described in Appendix [A] 
Overall, we find that each line could be decomposed into ^3 
components. The resulting emission line maps showing the 
flux, FWHM and radial velocity distributions are described 
in Section \3. II 



2.3 WIYN/SparsePak observations 

Complementary observations of NGC 253 were obtained 



with the SparsePak instrument (Bershady et al. 2004) on 



the WIYN 3.5-m telescope. SparsePak is a "formatted field 
unit" similar in design to a traditional IFU, except that 
its 82 fibres are arranged in a sparsely packed grid, with 
a small, nearly-integral core ( |Bershady et al.||2004| ). Each 
fibre has a diameter of 500 \im, corresponding to 4'! 69 on 
the sky; the formatted field has approximate dimensions of 
72 x 71.3 arcsecs, including seven fibres located on the north 
and west side of the main array separated by ~25 /; which are 
designed to sample the sky background. SparsePak is con- 
nected to the Hydra bench-mounted echelle spectrograph 
which at that time used a T2KA 2048 x 2048 CCD detector. 

We observed NGC 253 at two locations, one centred 
on the nuclear region at a=0 h 47 m 32!7, 6=-2b°17'19f.'0 
(J2000) (hereafter SP1) and one on the northern disk at 
a=0 h 47 m 29!5, £=_25°16 / 19 / . / 9 (hereafter SP2) on the night 
of 21st November 2005 with exposure times of 3x500 and 
3x1000 sees, respectively. The position of the SparsePak 
footprints (together with the VIMOS footprints) are shown 
in Fig.[l] overlaid on the WFI Ha image. Using the 860@30.9 
grating with an angle of 50.9° (giving a spectral coverage of 
6150-7050 A and dispersion of 0.44 A pix -1 ), we were able 
to cover the nebular emission lines of Ha, [N n] AA6548, 6583, 
and [S n]AA6716, 6731. A number of bias frames, flat-fields 
and arc calibration exposures were also taken together with 
the science frames. Since the seven "sky fibres" of the array 
fell on the galaxy itself, we took offset sky observations with 
exposure times equal to those of the corresponding science 
frames to facilitate sky subtraction. We therefore treated the 
"sky fibres" as normal science fibres and processed them ac- 
cordingly. 

Basic reduction was achieved using the CCDPROC task 
within the NOAO iraf package. Instrument-specific reduc- 
tion was then achieved using the hydra tasks also within 
the NOAO package. Cosmic-rays were cleaned by use of a 
simple minmax rejection as the individual frames were com- 
bined using imcombine. The separate offset sky observa- 
tions were then scaled and subtracted from the correspond- 
ing summed science frames. The final datafile contained 82 
reduced, wavelength calibrated and sky-subtracted spectra. 

In order to determine an accurate measurement of the 
instrumental contribution to the spectrum broadening, we 
selected two high S/N and isolated sky lines in one of the 
wavelength calibrated sky frames and fitted them with single 
Gaussians in all 82 apertures. We find the average instru- 



mental width at Ha is 1.77=b0.17 A or 81±5 kms -1 corre- 
sponding to a spectral resolution of R^3700. 

Again, the emission line profiles of Ha, [Nil] and [S n] 
were fitted with multiple Gaussians as described in Ap- 
pendix [A] Overall, we find that each line could be decom- 
posed into ^3 components. The resulting emission line maps 
showing the flux, FWHM and radial velocity distributions 
are described in Section! 



3 RESULTS 

3.1 VIMOS emission line maps 

In this section we present and describe the VIMOS emis- 
sion line maps created from line profile fits of the data de- 
scribed above. As mentioned above, each emission line in 
each spectrum was fit with a mult i- Gaussian model, result- 
ing in measurements of the flux, FWHM and radial velocity 
for each component (see Appendix [AJ. It is very important 
when fitting multi-Gaussian models to line profiles to assign 
each component to a particular map in such a way as to 
limit the confusion that might arise during analysis of the 
results, such as discontinuous spatial regions arising from 
incorrect component assignments. We assigned the bright- 
est component of Ha to CI, and the faintest to C2 (or C3 
if three components were identified). However, position 3 
made more sense if, where two components were detected, 
the redder was assigned to CI and the bluer to C2, regard- 
less of their fluxes. For [Sn], as with Ha, we assigned the 
brightest component in each case to CI, except for position 
3 where we assigned the redder to CI. No evidence of a third 
component was found in [S II] . 

Fig. [4] shows a selection of observed Ha+[Nll] line pro- 
files and best-fitting Gaussian models from a number of re- 
gions within our VIMOS fields (located with the correspond- 
ing letters in Fig. |5|. This sample represents the variety of 
different profile shapes we find across the IFU fields, and 
demonstrates the high quality of the spectra and the accu- 
racy of the line-fitting. Below each plot are the correspond- 
ing residuals, n, calculated using the following formula: 



fit _ data 
t/i t/i 



(i) 



where <7i are the uncertainties on yf ata . 

As mentioned in Section [2. 2[ the spectra in Q3 of each 
VIMOS field were blurred due to a focussing problem, mean- 
ing we could not perform a Gaussian decomposition of the 
profile shape. We therefore simply integrated the fluxes over 
each emission line of interest (Ha, [Nn]A6583, [Sll]A6717 
and A6731). The continuum level was then subtracted from 
each to recover the total line fluxes. 

Fig. [5] shows the line flux maps for the three identified 
Ha line components over the five IFU fields. The CI map 
(left panel) also includes the integrated line flux measured in 



Q3, where the focussing problems (see Section 2.2.1) meant 



that the line profiles could not be decomposed accurately. 
Since the spaxels in positions 1 and 2 are of a different size to 
those in positions 3, 4 and 5, we have made sure to present 
the fluxes in arcsec -2 units. Comparing Fig. [2^ and [H] re- 
veals that the depth of our IFU Ha spectroscopy and the 
continuum subtracted WFI Ha imaging is approximately 
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Figure 4. Example Hck+[Nii] line profiles extracted from the locations shown in Fig. [5] chosen to represent the main types of profile 
shapes observed over the field. Observed data are shown by a solid black line, individual Gaussian fits by dashed red lines (including 
the straight-line continuum level fit), and the summed model profile in solid red. Flux units are 10 — 16 erg s — 1 cm -2 A -1 spaxel - 1 . 
Below each spectrum is the residual plot (in units of cr; see text for further information), (a) An example line fit from position 1 showing 
a redshifted broad (FWHM~320 kms -1 ) component; (b) a profile from position 2 showing a blueshifted broad (FWHM~390 kms -1 ) 
component; (c) example of a clear split-line profile with velocity separation of ~260 kms -1 from position 3; (d) example of a triple 
component profile from position 4; (e) a spectrum from position 5 showing a broad component in [Nil] and not Ha; (f) example of a 
spectrum from position 4 with strong Ha absorption. 
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Figure 5. Ha flux maps in the three line components. The CI map includes the integrated line fluxes measured in Q3, where the focussing 
problems meant that the line profiles could not be decomposed accurately. A scale bar is given in units of 10 — 16 erg s — 1 cm -2 arcsec - 2 . 
The spaxels from which the example line profiles shown in Fig. [4] were extracted are labelled with the corresponding letters. 
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Figure 6. Ha FWHM maps in the three Ha line components. Left: CI, centre: C2 and right: C3. A scale bar is given in units of kms -1 , 
corrected for instrumental broadening. The vertical bands towards the bottom-left of position 5 are almost certainly instrumental in 
origin, although their precise cause has eluded us. 



the same. Both are considerably deeper than the available 
HST Ra imaging. 

FWHM maps for the three identified Ha line compo- 
nents are shown in Fig. [6j and the radial velocity maps 
are shown in Fig. [7] Where a comparison is possible, the 
velocities we measure are consistent with those measured 
by |Schulz fe Wegner| p992) ( given their significantly lower 
spatial resolution and the fact that they only fitted single 
Gaussians). Figs. [8] and [9] show position- velocity plots ex- 
tracted from a number of vertical (minor axis; PA=140°) 
and horizontal (major axis; PA=50°) l"-wide pseudo-slits 
whose positions and widths are shown with dashed lines 



in Fig. [7] Hereafter, we refer to these as horizontals^ or 
vertical##, where ## represents the y or x axis offsets of 
the pseudo-slit in arcsecs, respectively. The horizontal65 and 
vertical24 plots are in good agreement with those presented 
in |Prada et~aT] (1998) . 

Maps of the electron density distribution (as derived 
from the flux ratio of the [S 11] AA6717,6731 doublet, assum- 



ing T e = 10 4 K) are shown in Fig.llOj Before fitting the [S 11] 
lines, the data cube was smoothed in the spatial directions 
by a 3x3 rolling summation to increase S/N. A third line 
component could not be detected in [S 11] in any of the five 
fields, so only the CI and C2 maps are shown. 
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Figure 7. Ha radial velocity maps in the three line components. Left: CI; centre: C2 and right: C3. A scale bar is given in units of 
kms -1 , relative to ^ S ys- The dashed lines on the left panel indicate the positions and widths of a number of pseudo-slits from which 
the position-velocity data shown in of Figs. [8] and [9] were extracted. The solid and dashed white lines on the central panel highlight the 
approximate cone morphology determined by eye from the diverging Ha velocity components. 



The forbidden/recombination line flux ratios of 
[Sn](A6717+A6731)/Ha and [N n]A6583/Ha can be used 
as indicators of the radiation field strength and therefore 



the dominant ionization mechanism (Veilleux & Osterbrock 
|1987||Dopita fe Sutherland! 1995| [Dopita et al.|2000| |2006 f 



[Sll]/Ha is particularly sensitive to shock ionization be- 
cause relatively cool, high-density regions form behind shock 



fronts which emit strongly in [Sll] thus producing an en- 


hancement in [Sii]/Ha (jDopita 


1997 


Oey et al. 


2000). 


Figs. 11 and |12| show maps of the [Sll]/Ha and []N 


[ii] /Ra 



flux ratios, respectively, for the identified line components. 
Markers on the scale bar indicate values of log[Sn]/Ha = 
—0.4 and log[Nn]/Ha = —0.2, representing fiducial thresh- 
olds for the two ratios above which non-photoionized emis- 



sion is thought to play a dominant role (Dopita & Suther- 
land||1995| |Kewley et aLl|2001| jCalzetti et al.||2004| ). In en- 
vironments such as this, the most likely excitation mecha- 
nism after photoionization is that of shocks. Our line ratio 
maps are in agreement with those presented by Matsub- 
ayashi et ah] fl2009| ) and |Sharp fe Bland-Hawthorn| fl2010| ), 
given that we decompose our line profile shapes, and have a 
higher spatial resolution. 

In the following sections we describe what these maps 
and plots show, and comment on a number of the salient 
features in the area covered by our observations. Discus- 
sions of the meaning of some of these findings is presented 
in Section [4] 



3.2 Line shape distribution and general 
observations 

We detect at least two Ha line components over the ma- 
jority of the fields observed. In some regions in positions 
1 and 2, 4 and 5, we detect a third component, which, in 
most cases is narrower and/or fainter than the first two. 



Through examination of both the FWHM and radial ve- 
locity maps together (Figs. [6] and the Ha line profile 
can be classified into two main shapes/configurations: (a) a 
narrow component superimposed on a broader, fainter com- 
ponent at approximately the same velocity (e.g. Fig. [4^ or 
b); (b) a split narrow line where the two components are 
of approximately the same width (e.g. Fig. Type (a) is 
mainly found in position 2 just south of the nucleus (where 
we see very broad C2 widths; Fig. [6] middle panel), and 
position 1 just below and to the right of the missing Q3 sec- 
tion. Here we see C2 line widths of >400 kms -1 . This type 
is also found towards the lower-right (south) of position 1, 
but here the broad component is also the brightest and has 
therefore been assigned to CI. Type (b) is found over most 
of the rest of the fields, indicating the widespread presence 
of dynamically expanding material. The red and blue com- 
ponents vary in strength relative to one another (see below), 
and vary in width. In places, the two dynamical components 
both become very broad, reaching >300 kms -1 . Compared 
to positions 3 and 4, the C2 component in position 5 remains 
surprisingly broad over a large fraction of the field, where 
CI remains reasonably narrow (i.e. a type (a) narrow+broad 
configuration). However, unlike the type (a) shape found in 
position 2, here the broad C2 is clearly kinematically asso- 
ciated with the outflow cone structure (see below). 

In position 1, the wedge-shaped blueshifted region with 
its apex pointing to the nucleus (see Fig. [7] left panel) is 
indicative and consistent with the known disk rotation kine- 
matics - although no obvious evidence for the redshifted 
counterpart is seen in position 2, even after re-scaling the 
colour bar. 

Much of both the CI- and C2-emitting gas has [N n]/Ha 
and [Sn]/Ha line ratios (Figs. 11 and 12) that fall above the 



fiducial non-photoionization threshold, suggesting a signif- 
icant contribution to the ionization from shocks. Overall, 
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Figure 8. Position- velocity plots extracted from 1" wide pseudo-slits aligned with the galaxy minor axis at various positions along the 
dataset. The plots are labelled "vertical##" , where ## represents the x-axis intercept of the pseudo-slit on the map of Fig. [7] On these 
graphs, distance along pseudo-slit=0 (i.e. x=0) corresponds to y-axis=0 in Fig. [7] The velocities are given relative to v sys (250 kms -1 ). 
The vertical grey dashed line indicates the transition between positions 1/2 and 3/4. 



CI exhibits higher ratios, particularly in [Sn]. However, we 
also see localised regions of much higher line ratios, most 
likely representing pockets/filaments of strongly shocked gas 
(log([Sn]/Ha)>0.6 and log([Nn]/Ha)>0.7). These are cor- 
related with fainter line fluxes, a relationship that we have 



seen in similar environments in other galaxies (Westmo- 
quette et al.||2007a| |2009a ). 



Since we are observing the Ha-emitting gas superim- 
posed onto the background disk of the galaxy, stellar Ha 
absorption is expected, and indeed identified in some spax- 
els (see Fig. [4j). It is therefore possible that the regions of 
enhanced [N n]/Ha and [S ll]/Ha line ratios seen in positions 
3 and 4 are simply the result of Ha absorption affecting the 
line ratios. Given the spectral resolution and S/N of our 
data, it is very difficult to quantify the level of Ha absorp- 
tion in each spectrum. In order to investigate whether the 
regions of enhanced line ratios are due to Ha absorption, 
we have measured and mapped out the equivalent width of 
the integrated Ha line flux, and the continuum flux level 
near Ha. These maps are shown in Fig. |13| No correlation is 
seen between the regions of enhanced line ratios and either 



of these maps, leading us to conclude that the effect of Ha 
absorption is only slight, and our measured line ratios are 
real. 



The [S ii]-derived electron densities (Fig. 10) are high 
in the nuclear region, peaking at ~2000 cm -3 , and remain 
>700 cm -3 in the whole central disk. By a radius of ^10" 
the electron densities have dropped dramatically to close to 
the low density limit (few 100 cm -3 ). Localised regions of 
increased density are seen, with values up to 500-700 cm -3 . 
The density morphology does not follow any of the other 
patterns seen in the Ha kinematics. 



3.3 Outflow cone 

As described in the introduction, deep X-ray and Ha imag- 
ing have revealed that the inner ~1 kpc of the starburst out- 



flow within NGC 253 forms a limb-brightened cone (Strick- 



land & Stevens 2000). This structure can also be clearly 



identified both in our IFU maps and WFI Ha imaging. In the 
following, we describe its properties as derived from these 
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Figure 9. Position-velocity plots extracted from l /x -wide pseudo-slits aligned with the galaxy major axis at various positions along the 
dataset. The plots are labelled "horizontal##" , where ## represents the y-axis intercept of the pseudo-slit on the map of Fig. [7] On 
these graphs, distance along pseudo-slit=0 (i.e. x=0) corresponds to a>axis=0 in Fig. [7] Again, the velocities are given relative to v sys . 



observations. A simple kinematic model of the outflow cone 
is presented and discussed in Section [XT] 

The outflow cone is clearly seen in the morphology of 
the diverging Ha velocity components in the VIMOS radial 
velocity maps of Fig. [7| where the two components trace out 
a triangular shape extending south-east from the nuclear re- 
gion along the minor axis. We have attempted to highlight 
this formation, which has an opening angle of ^60-65°, with 
white lines on the C2 map. That the Ha cone is contained 
within the minor-axis spurs of H2 emission seen in the NIC- 
MOS image of Sugai et al. (2003), implies that molecular H2 



is entrained by the outflow. 

The P-V diagrams derived from the vertically-oriented 
pseudo-slits (Fig. [8} quite clearly show diverging Ha veloc- 
ity components emerging from a chaotic central region at 
disk heights of >15 /; . The coherence of the Ha velocities in 
the outflow region demonstrates that the 10 4 K, Ha-emitting 
gas layers must be quite thin and dynamically stable. The 
horizontal59 and horizontal 1 (Fig. |9| show how the cone 
expands from ~10" in width at a minor axis distance of 10" 
from the nucleus, to >40" at a projected distance of 30" . 
The morphology of the FWHM maps also broadly follows 



the cone shape, although the line widths vary dramatically, 
with many small-scale regions of much broader emission. 
Generally, the line widths remain above ^200 km s _1 in both 
CI and C2. The very broad components suggest there may 
be material filling the cone in some regions. 



3.4 Kinematically distinct features in the flow 

The flux, FWHM, and radial velocity morphologies are sug- 
gestive of small-scale, kinematically distinct, shell- like or fil- 
amentary features embedded within, and aligned with the 
outflow. For example, position 3 contains what appears to 
be a filament of material that is redshifted in both CI and 
C2 compared to the surrounding material (this shows up in 
the verticall6 P-V plot, Fig. [5] at x=50-55"). From the ra- 
dial velocity map, this feature appears to extend along the 
wind outflow direction, radially away from the nucleus. An- 
other region of redshifted material appears in C3 in position 
4. This region is again clearly kinematically distinct from 
the surrounding gas (see the horizontal41 P-V plot, Fig. [9] 
at x=35-40"). 

These distinct filaments/blobs of gas may be material 
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Figure 11. [S ll](A6717+A6731)/Hai flux ratio maps for CI (left) and C2 (right). C3 was not identified in [S n] . The marker in the 
colour-bar at —0.4 represents the fiducial ratio above which the excitation is likely to be dominated by non-photoionization processes 
(see text). 



in the foreground that is moving separately to the gas in 
the outflow. Alternatively, the redward velocity shifts may 
suggest that we are seeing through gaps in the foreground 
material, allowing the line-of-sight to penetrate either to a 
point deeper within the outflow, or gas associated with the 
background disk of the galaxy. 

These features are reminiscent of other systems with 



"bubbling winds" such as the LIRG NGC 3256 (Li'pari et al 



2004a see also Lfpari et al. 2004b). These type of flows 



seem to produce fragmented shells where the cool gas is al- 
ready accelerated, and thus differ somewhat from the radial 



streamers predicted to come from stationary clumps in a 
steady outflow (e.g. jPittard et al.|2005||Cooper et al.|2009 ). 



3.5 SparsePak emission line maps 

Figs. [14] and [15] show the SparsePak FWHM and ra- 
dial velocity maps, respectively. The [N n]A6583/Ha and 
[S n](A6717+A6731)/Ha line ratio maps are shown in 
Fig.[l6] As mentioned above, the footprints of the SparsePak 
observations are shown in Fig. [I] 

Two Ha line components are detected in the central 
disk and southern outflow regions. In the spaxels covering 
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Figure 12. [N n] A6583/Hai flux ratio maps in the three line components. Left: CI, centre: C2; right: C3. The marker in the colour-bar 
at —0.2 represents the fiducial ratio above which the excitation is likely to be dominated by non-photoionization processes (see text). 
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Figure 13. (a) Equivalent width of the integrated Ha line flux; (b) continuum flux level measured near Ha in units of 10 16 
erg s _1 cm -2 arcsec -2 . Some residual instrumental artefacts are evident around the field edges. 



the nuclear regions and a few associated with the south- 
ern outflow, the line shape is of the type (a) configura- 
tion (narrow+broad), with broad component widths of 200- 
450 kms -1 (Fig. 14). The remaining spaxels show type (b) 
line shapes (two split line components of ~same width). Af- 
ter running our tests and filters to remove incorrect or un- 
physical line fits, we assigned the brightest component to 
CI in all cases. A third line component can be identified in 
a few spaxels in the nuclear region. Here we assigned the 
additional narrower component to C3. 



As modelled by Sorai et al. ( 2000 ) and Das et al.| ( 2001 ) , 



the inner X2 orbits of the stellar bar extend diagonally over 
the full extent of position SP1 (as shown in Fig.[T]), meaning 



that the major axis dynamics are dominated by the inner 
bar region. 

Densities (maps not shown) of >300 cm -3 are found 
in places along the major axis (disk and starburst), but 
otherwise they remain close to or below the low density 
limit (<100 cm -3 ; as expected in the halo). Where the 
SparsePak and VIMOS observations overlap, the electron 
densities measured with both datasets are consistent (once 
the coarser spatial resolution of SparsePak is taken into ac- 
count). 



[Nil] /Ha and [Sn]/Ha line ratios (Fig. 16) are again 
consistent with those measured from our VIMOS data, 
where the coverage of the two datasets overlap. The highest 
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Figure 14. SparsePak Ha FWHM maps in the three Ha line components. Left: CI, centre: C2 and right: C3. A scale bar is given in 
units of kms -1 , corrected for instrumental broadening. 
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Figure 15. SparsePak Ha radial velocity maps in the three line components. Left: CI; centre: C2 and right: C3. A scale bar is given in 
units of kms -1 , relative to v sys - The VIMOS IFU footprints are overlaid. The diagonal dashed lines indicate the position and width of 
a pseudo-slit from which the position- velocity data shown in of Fig. |17| were extracted. 
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Figure 16. Left: SparsePak [N n] A6583/Ha flux ratio maps from summed line fluxes over all identified components; right: SparsePak 
[S ll](A6717+A6731)/Hoi flux ratio maps from summed line fluxes over all identified components. The marker in the colour-bars represents 
the respective fiducial ratios above which the excitation is likely to be dominated by non-photoionization processes. 
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Figure 17. Minor axis position- velocity plot including data from 
both the VIMOS and SparsePak datasets. Crosses represent ra- 
dial velocities extracted from the vertical24 VIMOS pseudo-slit 
shown in Fig. [8] whereas filled circles represent SparsePak veloci- 
ties extracted from a 10 /; -wide pseudo-slit aligned with minor axis 
as shown on Fig. ^] The offsets are calibrated such that x = 
corresponds to the nucleus. Velocities are corrected to the helio- 
centric reference frame, and given relative to v S y S (=250 kms" 1 ). 



line ratios are seen across the south-east of position SP1, 
coincident with the southern outflow cone. Line ratios fall 
dramatically just to the north-west of the nucleus, along a 
band parallel to the galaxy major axis. They then gradually 
rise again towards the north-west of position SP2, out to 
MOO" (>2 kpc) from the nucleus. 

To check consistency between our VIMOS and 
SparsePak datasets, we extracted the velocities along a 10"- 
wide pseudo-slit aligned with minor axis as shown on Fig. |15| 
and plotted them together with the VIMOS radial veloc- 
ities from the vertical24 pseudo-slit in Fig. |17| This plot 
shows how the measurements from the two datasets corre- 
spond very well. These velocities are also consistent with 
those measured by Schulz & Wegner (1992). 



4 DISCUSSION 

In this section we discuss the implications of our results on 
our understanding of the NGC 253 superwind. We then dis- 
cuss the relationships between the known properties of the 
central starburst, the disk, the molecular gas distribution, 
and the larger-scale properties of the aforementioned diffuse 
halo, to the observed characteristics of the inner wind. A 
number of the features that we refer to in the following sec- 
tions are illustrated schematically on Fig. [20] introduced at 
the end of Section |4.2[ to aid the reader in following our 
discussion. 



4.1 Kinematic modelling of inner wind cone 

We modelled the inner wind cone as described in Sec- 
tion |3.3| using the 3D morpho-kinematical modelling tool 
shape (Steffen et al. 2010). We restricted our fit to the re- 
gion between R= 15-35 from the nucleus (or ?/=35-55" on 



Fig. |7| where a coherent signature of an expanding conical 
outflow is observed. This shape therefore forms a conic frus- 
tum. Beyond this outer radius, the gap in our IFU coverage 
restricts our analysis, and it is not clear what happens to 
the cone structure here. 

We chose to model the particle velocities as flowing 
tangentially to the cone surface since this provided a good 
match to the observations, and was most easily understood 
physically in terms of material entrained by the high- velocity 
wind. Furthermore, a number of constraints to the model 
can be made from the observations. The fact that the two 
velocity components continue to diverge with radius implies 
that the outflow is accelerating. The increasing velocity of 
the redshifted component indicates that the far-side of the 
cone is angled away from the line-of-sight (assuming the gas 
flows tangentially to the cone surface). A smaller opening 
angle not only does not fit the morphology of the velocity 
maps of Fig. [8] but requires excessively large accelerations 
to fit the data. It was not required to change the inclina- 
tion of the cone from that of normal to the galaxy disk (i.e. 
i = 12°). 

After performing simulations over a range of frustum 
sizes, inclination angles, opening angles, and velocity laws, 
trying to match both the morphology of the Ha radial ve- 
locity maps, and the position- velocity diagrams of Fig. [8] 
we derived a best-fitting model whose parameters are given 
in Table [2] Fig. [l8| shows the model, both as seen projected 
onto the plane of the sky and side-on, where the arrows 
show the direction of the velocity vectors and illustrate how 
the flow speed increases with radius. Fig. [19] shows three 
of the observed P-V diagrams from Fig. [8] together with 
the simulated data from the best-fitting model, illustrating 
how well this simple model fits the data. There are indi- 
cations, however, that there are distortions from a simple 
cone shape and/or simple velocity law. For example, an off- 
set of —30 kms -1 needed to be applied to the model data in 
the vertical24 plot in order to match the observations. Addi- 
tionally, the velocity ellipse in horizontal41 (Fig.[9| does not 
appear to close on the left side, again indicating deviations 
from a simple geometry. 

The inferred diameter of the base of the outflow model 
(the top end of the cone) is ^20" (^380 pc). This, however, 
does not reflect the diameter of the energy injection zone 
(the starburst region) , but simply the width of the optically 
visible cone as it emerges from the veil of obscuration sur- 
rounding the nuclear regions. The best-fitting velocity law 
tells us that the Ha-emitting gas, as it enters this top end, 
is travelling at a deprojected outflow speed of 110 kms -1 , 
which then accelerates at a rate of 5z (where z is the non- 
projected physical distance along the cone in arcsecs). If the 
gas accelerates at a constant rate from its launch site, then 
the flow would become supersonic with respect to gas at 
10 4 K (sonic speed ^10 kms -1 ) almost immediately. 

By the time the gas reaches the end of the cone model 
(35"), it is travelling at 210 kms -1 . If this velocity law 
continues to hold to larger radii, it would be travelling at 
360 kms -1 at 50" (~1 kpc). This is much slower than the 
deprojected outflow speed of 390 kms -1 estimated for ma- 
terial at R < 16" by |Schulz fe Wegner| ( |1992| , a nd the 
~600 kms -1 outflow speed of the Ha gas in M82 (Shop- 



bell fe Bland-Hawthorn||1998). Since the escape velocity of 
NGC 253 is also estimated to be a few hundred kms -1 
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Figure 18. Best fitting kinematic model of outflow cone from 
Shape as seen (a) projected onto the plane of the sky, and (b) 
side on. The arrows in the top panel illustrate the direction and 
magnitude of the velocity vectors: along the cone walls, increasing 
with radius. 



(Heesen et al. 2009a| , it is not clear whether the Ka gas 
would be able to escape the galaxy's gravitational field. The 
hot X-ray emitting gas, however, must have an outflow speed 
which is faster than this, implying that the hottest material 
has the potential of escaping the galaxy. 

The best-fitting model has an opening angle of a — 60° 
between radii of 15-35", constrained quite tightly by the 
morphology of the Ha C2 velocity map (Fig. |7f. This is 
much larger than the 20-25° or 32° opening angles measured 
by |Strickland fe Stevens] ( |20QQ| ) and |Pietsch et al.|d2000| , re- 
spectively, although these were measured at larger radii (30- 
55"). This is also much larger than the ~25° opening angle 
of the M82 outflow ( |Shopbell fe Bland-Hawthorn|1998) . 

To measure the kinetic energy in the inner wind from 
our model, we assume the volume of the model conical 
frustum as described above is 2.1 x 10 7 pc 3 (assuming a 
cone wall thickness of 2" ~ 40 pc), and the electron den- 
sity in the cone walls is <100 cm -3 . This gives a mass of 
Ha-emitting gas of ^5 x 10 7 / M (where / is the fill- 
ing factor). To check this we can also apply the method 



X + ^ 
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Figure 19. Position- velocity plots from Fig. [8] (black +) with the 
best-fit model points overlaid (red x). The x-axis scale is the same 
as in Fig. [8] The vertical dashed line indicate the region within 
which we modelled the outflow cone. An offset of —30 kms -1 has 
been applied to the model data in the vertical24 plot in order to 
better match the data. 



of Pottasch (1980): an observed Ha surface brightness of 
~1 x 10 _ erg s _1 cm -2 arcsec -2 and electron density of 
100 cm - , results in a mass of 6xl0 6 M , which is m reason- 
able agreement if the filling factor /~0.1. Thus, assuming a 
mass of 10 7 M© and an outflow velocity of ^200 kms -1 , the 
Ha-emitting phase has a kinetic energy of ^4 x 10 54 ergs. 

The SNe rate in NGC 253 has b een variously esti- 
mated to lie between 0.03-0.3 yr" 1 (|Mattila & Meikle 



2001 



Ulvestad & Antonucci 1997). Following the argu- 



ments of Strickland et al. (2002), assuming the lowest SNe 



rate estimate gives a mechanical energy injection rate of 
^5 x 10 41 ?7therm erg s _1 (where rytherm is the SN energy 
thermalisation efficiency). Thus a starburst of even 5 Myr 
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Table 2. Parameters of the conic frustum model that best re- 
produces the observed kinematics of the southern outflow in 
NGC 253. 



Parameter 


Best-fit value 


Projected distance of frustum start, R\ a 


15" (280 pc) 


Projected distance of frustum end, R2 


>35" (>660 pc) b 


Diameter of the base of the outflow, D\ 


20" (380 pc) 


Opening angle, a 


60° 


Inclination angle, i 


12° 


Position angle (P.A.) 


140° 


Velocity law, v(z) a (kms -1 ) (") 


110 + 5z 



a Here we define R as the projected distance below the plane 
of the disk, and z as the physical distance below the plane in 
the ^-direction. 

b This distance is constrained by the limit of our IFU cov- 
erage. The cone is known to extend further out (see text for 
discussion). 



duration is able to produce enough energy (^ 8 x 10 54 ergs, 
with a conservative choice of ?7therm=10%) to power the ob- 
served Ha outflow in the inner wind. This result is in agree- 



ment with the findin gs of |Strickland &; Stevens (2000) and 
Matsubayashi et al.| ( |2009| . 



4.2 The origins and evolution of the inner wind 

The wind cone modelled above does not seem to be centred 
on the known nuclear starburst (located in the missing Q3 
section of position 2), but instead on a region to the north- 
east of the nucleus. This displacement is also evident in both 
the deep Ha image of iStrickland & Stevens (2000) and the 



NICMOS H 2 image of Sugai et al. (2003). What, then, is 



driving the outflow? If it is the nuclear starburst, how does 
the gas escape the nuclear region? What is collimating the 
flow? What are the conditions in the gas at the base of the 
flow? 

We can see from the vertically-aligned P-V diagrams of 
Fig. [8]that, in Ha, the outflow emerges at ~15" from a dy- 
namically chaotic nuclear region. That the extinction in the 
nuclear region is know to be high {Ay > 10; Kornei & Mc- 



Crady 2009 Fernande z- Ontiveros et al.|2009 ) indicates that 
within ~15 ' we see, in the optical, only the surface layers 
of the gas. This is consistent with the fact that we measure 
peak [S n] -derived electron densities of only ^2000 cm -3 , 
when observations of the high-densit y molecular gas tracer , 
HCN, imply densities of >10 5 cm" 3 ( |Knudsen et al.||2007| ). 
Drawing any conclusions about the gas conditions within the 
nucleus from optical observations is therefore not possible. 

IR and mm-wave imaging, however, shows that the 
molecular gas in the starburst region is distributed in a 
toroidal shape (Engelbracht et al.|1998||Garcia-Burillo et aT 
2000 Pa glione et al.|1995 2004 ), and kinematic analysis sug- 
gests that molecular gas is channelled inward through the 
action of the bar ( |Sorai et al.|2000 ) . This is presumably how 
the starburst is being fuelled, but how does it affect the out- 
flowing gas? The first piece of evidence comes from the SiO 
(v=0, J=2—>1) observations made with the IRAM interfer- 
ometer ( |Garcfa-Burillo et al.|2000| , and shown as contours 
overlaid on the HST F814W (/-band) image in Fig. [2}x The 
SiO emission reveals two nested molecular tori in the cir- 



cumnuclear disk that are aligned with the PA of the stellar 
bar. The location and size-scales of these tori suggest that 
they may play a significant role in confining or collimating 
the outflow on the smallest scales. In support of this idea, 



Tenorio-Tagle &; Munoz-Tunon (1998) discuss how the in- 
flow of matter to fuel the starburst can lead to a stationary 
hydrodynamic flow. In their model, the infalling material 
produces a ram pressure that easily exceeds that of the disk 
thermal pressure, and, when confronted with the mechanical 
luminosity from the starburst (i.e. the wind ram pressure), 
acts to channel the expanding starburst wind and collimate 
it into a kpc-scale, filamentary, biconical outflow. 

The second piece of evidence comes from the gas mo- 
tions seen in Br7 in near-IR VLT/SINFONI IFU observa- 



tions (Muller-Sanchez et al. in prep.; see also Miiller- Sanchez 
|et al.|[2010| ). These observations reveal that there are mul- 
tiple distinct outflow channels associated with the differ- 
ent components of the nuclear starburst, and that the out- 
flows first exit the nuclear region in an easterly direction, 
then quickly divert towards the direction of the minor axis 
(Muller-Sanchez et al. in prep.). This diversion or curving of 
the flow, also visible in our deep WFI Ha imaging (Fig.[2^), 
is supported by the changing orientation of the nuclear star- 
burst clumps and extraplanar dust lanes seen on the HST 
F814W image with distance. Within ^550 pc of the galaxy 
mid-plane, the dust lanes are oriented in an easterly direc- 
tion (particularly on the north-eastern side of the nucleus). 
Beyond this radius, the orientation of the dust lanes change 
to being consistent with the minor axis. Since the toroidal 
circumnuclear molecular disk(s) and the larger-scale stel- 
lar bar is oriented at PA « 70°, we hypothesise that these 
components have a strong influence on the collimation and 
direction of the flow in these inner regions, and that the hot 
wind flows not only down the pressure gradient but down 
any local gas density gradients as well. Given that the bar 
is rotating, the flow from within the nucleus must therefore 
change direction every few Myr (analogous to a garden sprin- 
kler), and it is not clear how this would allow a kpc-scale 
outflow structure to survive. 

A schematic overview of the salient points of NGC 253 
as discussed above is sho wn in Fig. |20| Her e we use the 
2MASS i^s-band image (Jarrett et al. 2003) as the basis 



to illustrate the main large-scale features of the galaxy, as 
this clearly shows the stellar disk, spiral arms and bar. Fur- 
ther to this, we include an inset, zoomed-in cartoon depict- 
ing the configuration of the nuclear regions and inner wind. 
This shows the circumnuclear molecular tori, and the on- 
sky orientation and configuration of the inner wind outflow 
as implied by the VIMOS Ha, SINFONI Br7 spectroscopy 
and the dust lanes seen on the HST F475W and F814W 
imaging. 

Finally, we know that the the superwind must entrain 
substantial amounts of cooler material (a certain proportion 
of which also mass-loads the hot flow) through the detection 
of PAH emission ( Tacconi-Garman et al.||2005 ) and far-IR 



filamentary dust emission extending out to 6-9 kpc (Kaneda 
|et al.|2009| . Furthermore, observations of near-IR H2 emis- 
sion with HST reveal sharply delineated filamentary struc- 
tures extending out to ^130 pc from the disk that follow the 



edges of soft X-ray emission (Sugai et al. 2003). This spa- 



tial coincidence suggests that the molecular emission orig- 
inates in an interaction zone between superwind fluid and 
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Figure 20. Labelled 2MASS Kg-band image of NGC 253 jJarrett et al.|2003| showing the main features of the galaxy relevant to our 
discussion and conclusions. The inset shows a zoomed-in cartoon depicting the configuration of the nuclear regions and inner wind. 



surrounding colder molecular gas, again suggesting entrap- 
ment of cool material. 



4.3 Evidence for a partial break-out of the cone 

The vertical33 P-V cut (Fig. |8| indicates that the two main 
Ha line components converge again at a distance of ^55" 
from the nucleus (x^20 // on Fig. [8}. However, the radial ve- 
locity map (Fig. [7) in the region where the vertical33 pseudo- 
slit cuts through position 5 suggests this convergence may 
only be a localised phenomenon. The point at which we see 
the convergence is at the limit of the limb-brightened cone 
region identified from the deep Ha images and soft Chan- 
dra X-ray emission by |Strickland fe Stevens] ( |2000| ) (Fig. [5]). 
This therefore indicates that we may be seeing a remnant 
shell of the outflow. 

However, the outflow is known to extend to >10' 
(>12 kpc) from the large-scale soft X-ray emission morphol- 



2002 ), implying that the wind must have flowed out past this 



inner region in the past. This suggests that the wind is only 
quasi-steady, and that successive waves of outflowing gas 
driven by distinct episodes of nuclear star formation (per- 
haps also linked to nuclear gas morphology and/or infalling 
gas kinematics) result in nested superbubbles that take time 
to rupture. At 55" our model predicts a deprojected outflow 
speed of 385 kins" 1 , giving a dynamical age of ^2.5 Myr, 
meaning that this bubble must have been driven by the most 
recent burst. Given the small number of massive star clusters 
that can be identified (~10-12), and their compact arrange- 
ment within the ^200 pc nuclear region, it would be quite 
surprising if the starburst did not have quiescent and active 
periods. 

Further evidence for the remnants of a ruptured su- 
perbubble may lie in the series of Ha clumps identified by 



ogy seen by ROSAT ( |Pietsch et aL]|2000| |Strickland et al. 



Strickland & Stevens (2000) in their deep Ha image (that 
they refer to as "additional clouds"), Matsubayashi et al. 
(2009] that they collectively refer to as Hn region 'B') and 



VIMOS-IFU observations of NGC 253 19 



our reduction of the ESO WFI Ha image (Fig. [2ji). These 
also clearly show up in the position 5 Ha CI flux map (al- 
though some of the knots fall in the missing Q3 section) . The 
corresponding CI velocity map shows no evidence for any ve- 
locity variations in these regions, meaning that the clumps 
are kinematically indistinct from their surroundings. (We do 
see spatially coincident redshifted gas in the C2 map, but 
this component is very faint in this region.) Mat subayashi| 
et al.| ( |2009) argue that these clumps (their Hn region B) 
are not part of the wind, but H n regions within the disk, be- 
cause their [Nil] /Ha and [Sll]/Ha line ratios are consistent 
with similar Hn regions in other disk galaxies. Further in- 
vestigation at higher spatial resolution and S/N is required 
to distinguish whether these clumps are really associated 
with the wind or simply part of the galaxy disk. 

4.4 The northern disk-halo interface 



diffuse UV, X-ray and far-IR emission are observed to reach 
out to a much larger scale (>8 kpc) perpendicular to the 



The SparsePak observations (Section 3.5) show that the Ha 



emission from the region to the north-west of the nucleus 
(north-western part of position SP1 and all of SP2) remains 
quiescent and well-behaved. Interestingly, no evidence for 
a northern counterpart to the southern outflow is found in 



Ha: only one component can be detected (Figs. 15 and 17), 
and its line widt hs r emain <80 kms -1 , with most below 
<50 kms -1 (Fig. 14). Radial velocities are instead consis- 
tent with simple disk rotation, with a smooth gradient being 
seen across the field between —40 and +25 kms -1 from the 
north-east to south-west (although this gradient does not 
show up on the map due to the colour bar scaling). Since 
there is clear evidence for the northern outflow in soft X- 
rays ( |Pietsch et al.|2000 ), why might there be no signatures 
of it in Ha? In the X-ray map, absorption by the inclined 
disk (the inclination of NGC 253 is such that the northern 
outflow is tipped away from the observer and the disk is 
in the foreground) dramatically reduces the X-ray emission 
flux up to R~3 kpc ( |Pietsch et al.p OQQ) from the nucleus. 
Since our SparsePak observations are well within this 3 kpc, 
we attribute the lack of any outflow signatures in Ha also 
to obscuration by the foreground disk. 

The gradual rise in line ratios towards the north-west 
of position SP2, out to >100 // (>2 kpc) from the nucleus 
(Fig. 16) is interesting. The rise is particularly prominent 



in [Sn]/Ha, although both ratios are high enough to be 
consistent with a significant contribution from shock excita- 
tion. However, since there is no kinematic evidence for the 
northern outflow, shock excitation would seem unlikely. It is 
known that diffuse ionized gas (DIG) in normal spirals can 
have high [Nil] /Ha and [Sll]/Ha ratios ( |Collins fc Rand| 
20011 lOtte et al.||2002|). That this is not due to shock heat- 



ing comes from high resolution studies of the Galactic DIG, 
which has the the same emission line ratio characteristics, 
but where very precise line widths can be measured (e.g. 
Reynolds et al. 1977 Reynolds 1985). Thus, we conclude 
that the high line ratios result from excited DIG in the halo, 
with some contribution from underlying disk emission. 

4.5 The connection between the inner wind and 
the kpc-scale X-ray superbubbles 

As mentioned previously, apart from the nuclear outflow (in- 
ner wind) focussed on in this paper, "horn-shaped" lobes of 



disk (|Pietsch et al.|2000l [Strickland et al.|2002| |Bauer et al 
2008; |Kaneda et al.||2009| . In this section we briefly discuss 
how the two components of the outflow may be connected. 
Does the nuclear starburst-driven inner wind flow examined 
in this paper power the large-scale lobes, or are they formed 
through another mechanism? 

Assessing all the evidence, it is not clear to us that there 
is a link between the obviously nuclear starburst-driven in- 
ner wind and the large-scale X-ray lobes. Firstly there ap- 
pears to be a large disparity in scales between the colli- 
mated, conical inner wind flow and the widely spaced X-ray 
"horns" , suggesting an unphysically rapid expansion in the 
radial direction with increasing vertical height. Without fur- 
ther confinement, it is not clear how an X-shaped morphol- 
ogy could form. Recent observations of the synchrotron halo 
of NGC 253 suggest that a cosmic-ray (CR) driven "disk- 
wind" may play a significant part in driving material out 



into the halo and create the large-scale lobes (Heesen et al 



2009a b). The fact that these authors find a remarkably con- 
stant vertical CR bulk speed over the entire extent of the 
disk shows that the CRs do not originate from a localised 
(centrally located) source, but from the whole disk. Further- 
more, the vertical bulk speed of the CR-driven gas is large 
(^280 kms -1 ), meaning the kinetic energy of CR-driven gas 
is certainly high enough to expel material out of the galactic 
plane. 

On large-scales, the magnetic field structure of NGC 253 
resembles the X-shape of the Ha and soft X-ray lobes 
([ Strickland et al.||2002| |Heesen et aLl|2009a| ). |Heesen et aT 



2009b| found that both a disk wind model or a superwind 
model in conjunction with a large-scale dynamo action could 
explain the X-shaped halo magnetic field. However, they 
argue, such X-shaped magnetic field structures have been 
observed in several edge-on galaxies, not all with a central 
starburst and superwind. This suggests that the X-shaped 
halo field might be due solely to the disk wind. The super- 
wind may then be (partly) collimated by the halo field on 
large-scales, or perhaps the superwind may act further to 
amplify and align the halo field by compression of the lobes 
of the expanding superbubbles. 

As mentioned in the introduction, Boomsma et ah] 
( |2005| ) interpret the presence of a substantial mass of ex- 
traplanar Hi as a result of stellar feedback. Hee sen et al.| 
(2009b) argue that the asymmetric distribution of the gas 



makes a minor merger origin very unlikely, and go on to 
surmise that a disk wind is a good candidate for lifting this 
Hi gas out into the halo. However, to lift that amount of 
material to such large distances implies that a significant 
injection of energy would be needed (~5 x 10 55 erg) over 
a reasonably long time (~50 Myr). The disk does indeed 
contain a great many compact, bright Hn regions as shown 
in our WFI Ha image, and the existence of strong radio 



emission outside the 200 pc starburst region (Carilli et al 
1992) |Beck et aL||l994| |Carilli|[l996| |Heesen et aH ~ 



suggests a significant amount of star formation. Ulvestac 



2009a) 



(2000 1 estimates that the region outside the central star- 
burst may account for at least 20% of the recent star forma- 
tion in NGC 253. Energetically, therefore, it seems feasible 
that the plumes are powered by a disk wind and not the 
starburst-driven outflow. Recent, deep imaging by [Davidge 
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(2010) show that the north east portion of the NGC 253 



disk seems to have experienced excess star formation for the 
last few lOs-lOOs Myr. Although the 'disk-wind' mechanism 
was not mentioned specifically, the suggestion of a causal re- 
lationship between the elevated SFRs and the gas emission 
through feedback from a 'boiling disk' (Sofue et al. 1994) 
was made. 



5 SUMMARY 

In this paper we have presented high resolution optical in- 
tegral field unit (IFU) observations and associated deep Ha 
imaging of ionized gas in one of the nearest nuclear starburst 
galaxies with a large-scale superwind, NGC 253 (3.9 Mpc), 
in an attempt to understand the outflow in more detail. Our 
main results and conclusions can be summarised as follows: 

• We confirm the presence of an expanding outflow cone 
between radii of 280-660 pc from the nucleus, oriented along 
the minor axis. This structure is well-defined in the Ha mor- 
phology and kinematics. Dynamical modelling indicates a 
wide opening angle (~60°), an inclination consistent with 
being normal to the disk (12°), and outflow speeds of a few 
100 krns" 1 . 

• The [Nil] /Ha and [Sll]/Ha line ratio maps suggest a 
great deal of the gas in the wind is shocked, with localised 
regions of high line ratios most likely representing pock- 
ets/filaments of strongly shocked gas. Ha line widths of 
>400 krns" 1 (FWHM) are measured in regions associated 
with the wind cone, suggesting that there is material filling 
the cone in some regions. 

• By comparison to VLT/SINFONI IFU near-IR emis- 
sion line data (Muller-Sanchez et al. in prep.), we can trace 
the outflow back to the nuclear driving sources, and find ev- 
idence for multiple distinct outflow channels implying that 
the superwind is driven by more than one starburst region. 
From the kinematics, the cone appears partially closed in 
one place. These two pieces of evidence together strongly 
suggest that the wind is only quasi-steady, and that succes- 
sive waves of outflowing gas have been driven by distinct 
episodes of nuclear star formation to form the structure we 
see today. 

• The location and size-scales of the nested tori seen in 
SiO emission by Garcia- Bur illo et al. ([2000) suggest that 
they may play a significant role in initially confining or colli- 
mating the outflow. It is thought that molecular gas is chan- 



nelled inwards through the action of the bar (|Sorai et al. 



2000}. Thus, the ram pressure of this infalling material act- 
ing against the ram pressure of the wind outflow may also 
act to collimate the flow somewhat, as discussed by |Tenorio-| 
Tagle & Munoz-Tunonl (11998). 



• The SINFONI near-IR emission line measurements 
show that the outflow first exits the nuclear region in an east- 
erly direction, then quickly diverts towards the direction of 
the minor axis. This deflection, also visible in our deep WFI 
Ha imaging, is supported by the changing orientation of the 
extraplanar dust lanes seen on the HST F814W image with 
distance. Since the toroidal circumnuclear molecular disk(s) 
and the larger-scale stellar bar is oriented at PA « 70°, we 
hypothesise that these components have a strong influence 
on the direction of the flow in these inner regions, and that 



the hot wind flows not only down the pressure gradient but 
down any local gas density gradients as well. 

• No evidence is found in our optical Ha SparsePak ob- 
servations for an outflow on the north-western side out to 
>2 kpc. We surmise this is due to obscuration by the fore- 
ground disk (c.f. jPietsch et al.|2000 ). 

• The connection between the inner (R<1 kpc) Ha-bright 
conical outflow focussed on in this paper and the large-scale 
(R<10 kpc) X-ray "horns" is also discussed in the context 
of whether or not the inner starburst-driven wind powers 
the large-scale lobes. Heesen et al. ( 2009a|b ) find evidence 
for a cosmic-ray wind driven by the disk itself at speeds of 
^280 krns -1 . We conclude that the X-ray lobes are powered 
by a cosmic-ray driven disk wind (Heesen et al. 2009a|b ). 
Since the disk contains a great many compact, bright Hll 
regions as shown in our WFI Ha image, and strong radio 
emission from the disk suggests a significant amount of star 



formation is occurring outside the central s tarburst ([C arilli 
et al.||1992| |Beck et al.|^994| |Carilli|[l996] |Ulvestad||2000[ 
Hees en et al.||2009a| , we favour the explanation that the 
large-scale plumes are powered by a disk wind, with, if any- 
thing, only a small contribution from the starburst-driven 
inner wind. 
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APPENDIX A: DECOMPOSING THE 
EMISSION LINE PROFILE SHAPES 



Following the methodology first presented by Westmoq uette| 
et al.| ([2007a), we fitted multiple Gaussian profile models 
to each emission lin e see n in the VIMOS (Section 2.2) and 
SparsePak (Section 2.3) datasets using an idl x fitting 



package called pan (Peak ANalysis; Dimeo 2005) to quan- 
tify the gas properties observed in each IFU field. The high 
S/N and spectral resolution of these data have allowed us 
to quantify the line profile shapes of these lines to a high 
degree of accuracy. 

Each line in each of the 1200 spaxels (VIMOS IFU 
less Q3) and 82 spaxels (SparsePak IFU) was fitted using a 
single, double, and triple Gaussian component initial guess 
(no lines were found to need more than three components). 
Line fluxes were constrained to be positive and widths to be 
greater than the instrumental contribution to guard against 
spurious results. 

Although 15000 long slit echelle spectra taken by 
one of us (JSG) suggest some differences may exist between 
the Ha and [Nil] emission line profiles, these are not de- 
tected in our lower resolution data. Therefore, to fit the Ha 
and [Nn]A6583 lines, for each component we constrained the 
wavelength difference between the two Gaussian models to 
be equal to the laboratory difference, and the FWHMs equal 
to one another. This constrained approach, also taken for fit- 
ting the [S n] doublet, dramatically improves the quality of 
the fits, particularly for spectra with reduced S/N. Multi- 
component fits were run several times with different initial 
guess configurations (widths and wavelengths) in order to 
account for the varied profile shapes, and the one with the 
lowest x 2 ^ statistic was kept. However, we note that the 
X 2 minimisation routine employed by PAN is very robust 
with respect to the initial guess parameters. 

To determine how many Gaussian components best fit 
an observed profile (one, two or three), we used a likelihood 
ratio test to determine whether a fit of n+1 components was 
more appropriate to an n-component fit. This test says that 
if the ratio of the x? statistics for the two fits falls above 
a certain threshold, then the fits are considered statistically 
distinguishable, and the one with the lower % 2 can be se- 
lected. Here we determine the threshold ratio by visual in- 
spection of a range of spectra and fits. This method is a more 
generalised form of the formal F-test that we have used in 



previous work, e.g. Westmoquette et al. (2007a 2009a). We 
have chosen to modify our methodology in light of work we 
have recently become aware of that cautions against the use 
of the F-test to test for the presence of a line (or additional 
line component; |Protassov et al.|2002| ). In fact these authors 
show that using any kind of likelihood ratio test in this kind 
of situation is statistically incorrect, but given the absence of 
a statistically correct alternative that could be applied sen- 
sibly to the volume of data presented here, we have chosen 
to opt for this generalised x 2 ratio test approach. 

This test, however, only tells us which of the fits (sin- 
gle, double or triple component) is most appropriate for the 
corresponding line profile. Experience has taught us that we 
need to apply a number of additional, physically motivated 
tests to filter out well- fit but physically improbable results. 

Another addition to our methodology was to bootstrap 
the uncertainties on the fit parameters using Monte Carlo 



techniques, assuming that the noise is uncorrelated. Using 
the fit results determined from the above step as the initial 
guess, we re-fitted each profile 100 times, each time seeding 
the spectrum with randomly generated noise based on the 
standard deviation of the continuum. The standard devia- 
tion on each parameter (flux, FWHM and radial velocity) 
over the 100 results was then adopted as the associated un- 
certainty for that parameter, while the mean of those results 
was adopted as the actual value. Uncertainties range from 
1-30% depending on the S/N of the spectrum and how many 
components are fit. 

The final step of the fitting procedure was to assign each 
fit component to a particular map in such a way as to limit 
the confusion that might arise during analysis of the results, 
such as discontinuous spatial regions arising from incorrect 
component assignments. How we assigned the components 
in each case is described in the relevant results sections (3.1 



and 3.5 ) 



